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Abstract 



The decays D* — > Dtt and D* are well described by heavy meson 

chiral perturbation theory. With the recent measurement ol B{D*^ — > D^j), 
the D*^, D*^, and D* branching fractions can be used to extract the D*Dtt 
and D*Dj couplings g and /?. The D* D'j decays receive important cor- 
rections at order ^/m^ and, from the heavy quark magnetic moment, at order 
l/nic- Here all the decay rates are computed to one- loop, to first order in 
ruq and l/rUc, including the effect of heavy meson mass splittings, and the 

counterterms at order niq. A fit to the experimental data gives two possi- 
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ble solutions, g = 0.27l;g^ f3 = 0.851;? 1;? GeV"! or g = 0.76 + °3+-2 



p = A.90ti i^7° GeV^i . The first errors are experimental, while the second 
are estimates of the uncertainty induced by the counterterms. (The experi- 
mental limit r£)*+ < 0.13 MeV excludes the g = 0.76 solution.) Predictions 
for the D* and B* widths are given. 
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I. INTRODUCTION 



Combining chiral perturbation theory with heavy quark effective theory (HQET) gives a 
good description of the low energy strong interactions between the pseudo-goldstone bosons 
and mesons containing a single heavy quark. Due to heavy quark symmetry (HQS) there 
is one coupling, g, for 0*071, 0*0*71, B*B7i, and B*B*7i, and one coupling, j3, for 0*0j, 
0*0*'-f, B*B'y, and B*B*^ at leading ordeiQ. The value of the coupling g is important, 
since it appears in the expressions for many measurable quantities at low energy. These 
include the rate B 0^*^7i£i'i form factors for weak transitions between heavy and light 
pseudo-scalars decay constants for the heavy mesons weak transitions to vector 

mesons ||^, form factors for B 0^*Hi'i and heavy meson mass splittings (for a 
review see However, the value of g has remained somewhat elusive, with numbers in 

the literature from ~ 0.2 to 1.0. Recently, a CLEO measurement [l^ of O*^ has 
brought the experimental uncertainties to a level where a model independent extraction of 
g is possible from O* decays. 

As a consequence of HQS the mass splitting between O* and O mesons is small (of order 
Aqcj^/mc), leaving only a small amount of phase space for O* decays. In the dominant 
modes, O* —>■ O71, and O* O'j, the outgoing pion and photon are soft making the chi- 
ral expansion a valid framework. The branching ratios for 0*~^ decay are O'^ 7C^ (67.6%), 
L>+7rO (30.7%) and D+7 (1.7%) [g. A can only decay into D° 7r° (61.9%) and 
i5°7 (38.1%) [^] since there is not enough phase space for O^ . The O* decays pre- 
dominantly to Og'-f (94.2%) with a small amount going into the isospin violating mode 
OgTc^ (5.8%) |]T^. Since a measurement of the widths of the O* mesons has not yet been 
made, it is only possible to compare the ratios of branching fractions with theoretical pre- 
dictions. The ratio R+ = B{0*+ 0^7i+)/B{0*+ 0+7i°) is fixed by isospin to be 

= 2|A;^+p/|A;^op = 2.199 ± 0.064 |jl2[ (where fc^+,o are three momenta for the outgoing 
pions in the O* rest frame). This value is often used in experimental extractions of the 
branching ratios to reduce systematic errors. 

It is interesting to note that the quark model predictions |M for O*^ and O*^ decays 



^Where it is meaningful we use vr to denote any member of the pseudo-goldstone boson SU(3) 
octet, and D* and D for any member of the triplets (D*^ , D*^ , D*) and {D^, D^,Ds) with a similar 
notation for B* and B. 
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agree qualitatively with the data. One can understand, for instance, why the branching 
ratio B{D*^ D^l) is small compared to B{D*^ —>■ D^'-f). In the quark model the photon 
couples to the meson with a strength proportional to the sum of the magnetic moments of 
the two quarks, fi2 = 2/(3 r/ic) — 1/(3 m^) for D*^ and /ii = 2/(3 mc) +2/(3mu) for 

jj*o _^ £)0^. The rate for the former is then suppressed by a factor 



mu/nic + ly 

where we have used mass ratios appropriate for constituent quarks, mu/md ^ 1, md/rric — 
mu/rric ~ 1/4. This suppression results from the opposite signs in fii and /i2, which in turn 
follow from the (quark) charge assignments and spin wavefunctions for the heavy mesons. 

In the quark model g = 1 and /3 ~ 3 GeV~^, while for the chiral quark model g = 0.75 
15| . Relativistic quark models tend to give smaller values, g ~ 0.4 [|T6|, as do QCD sum 



rules, ^ ~ 0.2 - 0.4 1] 

Our purpose here is to use heavy meson chiral perturbation theory at one-loop to extract 
the couplings g and (3 from D* decays. In other words, we wish to examine how sensitive 
a model independent extraction of g and f3 is to higher order corrections. For D* D'y, 
analyses beyond leading order have included the heavy quark's magnetic moment which 
arises at l/m^ [|r^,^, and the leading non-analytic effects from chiral loops proportional to 
■^ymq y/rriq terms proportional to both mx and m^r were found to be important. These 
effects do not introduce any new unknown quantities into the calculation of the decay rates. 
For D* D'y and the isospin conserving D* Dn decays the effect of chiral logarithms, 



In (/i/m^), have also been considered pOf. These are formally enhanced over other 
corrections in the chiral limit, rriq — > 0, however, the choice of the scale fi leads to some 
ambiguity in their contribution. (This scale dependence is cancelled by unknown couplings 
which arise at order rriq in the chiral Lagrangian.) The isospin violating decay D^t^^ 



has only been considered at leading order, where it occurs through i] — tt mixing |21 



In this paper the investigation of all D* decays is extended to one-loop, including sym- 
metry breaking corrections to order rUg and l/rric. Further l/rric and rriq contributions 
considered here include the effect of nonzero D*-D and D^-D^ mass splittings, and the 
exact kinematics corresponding to nonzero outgoing pion or photon energy in the loop di- 
agrams. (Their inclusion is motivated numerically since m^o ~ uid* — m£, ~ m£,^ — mo, 
and the decay D* Dn^ only occurs if m^i* — m^) > m^o.) To simplify the organization of 
the calculation these splittings will be included as residual mass terms in our heavy meson 
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propagators. This gives new non-analytic contributions to the D* Dtt^ and D* 
decay rates. (To treat the mass sphttings as perturbations one can simply expand these 
non- analytic functions.) At order niq there are also analytic contributions due to new un- 
known couplings which are discussed. These new couplings can, in principle, be fixed using 
other observables. We estimate the efi^ect these unknown couplings have on the extraction 
of g and (3. 

The calculation of the decay rates to order and l/rric is taken up in section II. In 
section III we compare the theoretical partial rates with the data to extract the D*Dtt and 
D*D'y couplings and discuss the uncertainty involved. Predictions for the widths of the D* 
and B* mesons are also given. Conclusions can be found in section IV. 



II. DECAY RATES FOR D*^, L>*+, AND D* 

In this section we construct the effective chiral Lagrangian that describes the decays 
D* Dtt and D* —>■ D'-f to first order in the symmetry breaking parameters rrig and l/rric. 
The eight pseudo-goldstone bosons vr* that arise from the breaking SU{3)l x SU{3)j^ 
SU{3)v are identified with the pseudoscalar mesons {tt^ ,7t~^ ,7t~ ,1]) . These 
can be encoded in the exponential representation E = = exp(2i7r'AY/), where A* are 
3x3 matrices such that 





/ nyV2 + 7]/VQ 




K+ \ 




71 A = 




-7iyV2 + T]/V6 




(2) 













and f = 130 MeV. For the triplets of heavy mesons {D°, D+, and D*+, 

D^) we use the velocity dependent fields Pa{v) and -Pa^(f ) (a=l,2,3) of HQET. These are 
included in a 4 x 4 matrix which transforms simply under heavy quark symmetry 

Ha = ^^[P:'l,-Pa7,], (3) 

and satisfies '^Ha = Ha = —Ha'^. Including the quark mass term niq = diag(mu, rrid, rris) the 
lowest order Lagrangian is then 

Co = ^Tr d^T. d^T) + :^!^Tr(m,S + m,St) - Tr H^iv ■ D^^H, + gTi H^H.^^^.Al , (4) 
o 4 

where the derivative = Sabd^ — V^b, and Ha = l^H\'-j^. The vector and axial vector 
currents, Va^ = \{S}d^^ + ^d'^^'^) and A'^^ = ^{i'^d^^ — ^S'^^^), contain an even and odd 



number of pion fields respectively. The Lagrangian in Eq. (|) is invariant under heavy quark 
flavor and spin symmetry. It is also invariant under chiral SU (3) l x SU (3) r transformations, 
where E — LTiR\ ^ ^ L^U'^ = U^R\ H — > HU'^ ^ if we take the quark mass (which breaks 
the chiral symmetry) to transform as — > Lmq 

The last term in Eq. (Q) couples P*Pt{ and P*P*'k with strength g and determines the 
decay rate D* —>■ Dtt at lowest order. Going beyond leading order involves including loops 
with the pseudo-goldstone bosons, as well as higher order terms in the Lagrangian with more 
powers of rUg, l/rric, and derivatives. At order rrig ~ l/rric the following mass correction 
terms appear 

Cm = T^Tr Haa^'Haa^, + 2AiTr HaHbml + 2X[Ti HaHaml, (5) 

where = ^{^rrig^^ + C.'^mg^). The A'^^ term can be absorbed into the definition of niff 
by a phase redefinition of H. The A2 term is responsible for the D*-D mass splitting at 
this order, A = m£,* — rriu = ~2X2/mc. The term involving Ai splits the mass of the 
triplets of D and D* states. Ignoring isospin violation this splitting is characterized by 
6 = rriD* — mo* = frtOs ~ = 2Ai(ms — m) where rh = tUu = rrid. For the purpose of 
our power counting 6 ~ nig ~ 1/itLc ~ A. The effect of these mass splitting terms can 
be taken into account by including a residual mass term in each heavy meson propagator. 
Since we are interested in decay rates we choose the phase redefinition for our heavy fields 
to scale out the decaying particle's mass. For D*^ and D*^ decays the denominator of our 
propagators are: 2v ■ k for and 2{v ■ k — 5) for D*, 2{v ■ k + A) for and D+, 

and 2{v ■ k + A — 6) ioi Dg. For the D* decays the denominators are the above factors plus 
26. (If we scaled out a different mass then the calculation in the rest frame of the initial 
particle would involve a residual 'momentum' for the initial particle, but would yield the 
same results.) This results in additional non-analytic contributions from one- loop diagrams 
which are functions of the quantities A/m^. and Formally, m^. ~ rUg ~ A ~ 5 and 

one can expand these contributions to get back the result of treating the terms in Eq. (^ 
as perturbative mass insertions. 

Another type of l/rric corrections are those whose coefficients are fixed by velocity repa- 
rameterization invariance [^,0] 



SC, = -77^ Tr Ha{tD)lH, + -^i:iH,{iD^^^v ■ A,a - tv ■ A^c'^'^jH^'y^'y^ • (6) 
The first term here is the HQET kinetic operator, Okm = 2m^^^ (^-^)^ ^v, written in terms of 



the interpolating fields Pa and P*'^. In conjunction with the HQET chromomagnetic opera- 
tor, Omag = 2m^^^ ^ ^^apG^^ hy, these Contributions to the Lagrangian modify the dynamics 
of the heavy meson states. They give l/rric corrections in the form of time ordered products 
with the leading order current , which induce spin and flavor symmetry violating correc- 
tions to the form of the D*Dtt coupling. We account for these corrections by introducing 
the couplings gi and (72 in Eq- (0) below. The last term in Eq. (|^) contributes at higher 
order in our power counting since it is suppressed by both a derivative and a power of l/rric- 
Further terms that correct the Lagrangian in Eq. @) at order niq ~ 1/ uic include 00 



A2 ^^a^^u ifi iD^^tja''"cc ' a 2 ^^u^^ a i ^j, i ^bc^cb 

X X 

+ ^TiHaHb^^^.tv ■ DbcA>^,a + j-TiHaHb^^^^iDlv ■ A,, 

+ ^Tr HaHbi.i^Al + ^Tr H^i^i^HbAl + . . . , (7) 

ruQ 

where D^^A^^ = d^A^^ + [V"', A^]ba and = 47r/. The ellipses here denote terms linear in 



^ = ll^'^q^^ — ^^^gO which contribute to processes with more than one pion, as well as 
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terms with {iv ■ D) acting on an H. For processes with at most one pion and H on-shell the 
latter terms can be eliminated at this order, regardless of their chiral indices, by using the 
equations of motion for H. The Hi coefficients contain infinite and scale dependent pieces 
which cancel the corresponding contributions from the one-loop D* Dtt diagrams. For 
the Ki and n'l terms only the combination ki = Ki + k'^ will enter in an isospin conserving 
manner here. (The combination Ki — k,[ will contribute an isospin violating correction to 
i?^.) At a given scale fi, the finite part of ^3 can be absorbed into the definition of g. The 
decays D* Dtt have analytic contributions from ki and K5 at order m^. 

For rriQ = rric the term in Eq. (0) involving gi can be absorbed into g (this term only 
enters into a comparison with B* decays). The term g2 breaks the equality of the D*D7t and 
0*0*71 couplings. Since we only need the coupling D*D*tt in loops we can also absorb g2 
into the definition of g. Thus, our g is defined as the D*D7i coupling with l/mq corrections 
arising in relating it to the couplings for D*D*-k and B^*^ B*-k. 



^The k'i term was not present in The factor Bq/K^ is introduced here for later convenience. 



6 



a) b) c) 

/ \ / \ / \ 
I 1 I 1 I 1 

D D* D D* D D* D* D* D* 

FIG. 1. D and D* wavefunction renormalization graphs. The dashed hne represents a 
pseudo-goldstone boson. 

The terms in Eq. (|^) involving 62 and ^3 contribute to D* Dtt^, entering in a fixed 
linear combination with the tree level coupling g of the form g — {62 + 63)v ■ k/A^. These 
are ~ 10% corrections for the decays D* Dn. The energy of the outgoing pion is roughly 
the same for all three decays, f ■ A; ~ .144 GeV. Therefore, it is impossible to disentangle the 
contribution of 61^2 from that of g for these decays, and the extraction of g presented here 
will implicitly include their contribution. For other processes involving pions with different 
V ■ k these counterterms can give a different contribution. This should be kept in mind when 
this value of g is used in a different context. 

Techniques for one-loop calculations in heavy hadron chiral perturbation theory are well 
known and will not be discussed here. Dimensional regularization is used and the renormal- 
ized counterterms are defined by subtracting the pole terms 1/e — 7 + log (47r). The decays 
£)*o ^ £,0^0 g^^^ _^ 2^+^o^ g^^^ ^* _^ ^^^0 ^^^^ jgpj^y ^j^^gg rl, and Tl given by 



pa 9^ 



12 



7°- 



2 

I zTalS 



(8) 



Here fc" is the three momentum of the outgoing pion, Z'^ contains the vertex corrections, 
and = Zf), Z1) contains the wavefunction renormalization for the D* and D. When the 
ratio of to the D* —>■ D'-f rate is taken Z^^ will cancel out. However, Z^^ does contribute 
to our predictions for the D* widths, where the ratio Z^^/Z^ will be kept to order g^. The 
graphs in Fig. |l| give 

2 2 

^D = ^ + 7T^(Aj.,A;i)|[3m,2 - 6(A + do)'] log (^) + 3 G^{m„ A + do)] , 

2 2 

Z%.=1 + j£jY2iKbKi){[^ra^ - < - 2(rfo - A)^] log + 2 Gi(m„ do) 

+ GiK,rfo-A)j, (9) 



where is the mass of vr*, do = S^^6 for D*^ and D*^ decays and do = {S^'^ — 1)6 for D* 
decays. The notation in Eq. (P) assumes that we sum over 6 = 1,2,3 and i = 1, . . . , 8. The 
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FIG. 2. Nonzero one-loop vertex corrections for the decays D*^ D^n^ and D*^ D^tt^ 
(a,b,c) and the pseudo-goldstone boson wave function renormahzation graph (d). 



logarithms agree with |^ , except that we have kept terms of order ~ in the prefactor 
since these terms are enhanced for niq 0. Analytic terms of order ~ are neglected 
since they are higher order in our power counting. The function Gi{a, b) in Eq. (P) has mass 
dimension 2. It contains an analytic part proportional to a^, and a non-analytic part which 
is a function of the ratio b/a. The expression for Gi can be found in the Appendix. 

For a = 1, 2 the decay proceeds directly so that at tree level Z^f /Z^'^ = 1. At one loop 
we have non-zero vertex corrections from the graphs in Fig. Ha,b,c. As noted in pO[, the 



two one-loop graphs that contain a D^*^D*mT vertex (not shown) vanish, and the graph in 
Fig. lie cancels with the 7r° wavefunction renormahzation in Fig. |^d (this is also true for 
D*^ D^Ti^ and D* D^n^). Therefore for a = 1,2 the vertex corrections are 



„2 \i \1 \i\ ( ,,2 
9 ^ab\b\a I 1 / 



2 



log (^) + -{-dl + did2 + dl- 2di do - 2 dl) 



3 



+2Fi(mi,rfi,d2) -4Fi(mi,rfi,do)| + Ql{,hi,K^) , (10) 

where here d^ = 6^^6, di = k ■ v + do, d2 = —A + do, and k is the outgoing momentum 
of the 7r°. The coefficient of the log (yU^/mf ) term agrees with |2^. The function Fi 
has mass dimension 2 and contains both analytic and non-analytic parts, g'^^^ contains the 
dependence of the rate on the (renormalized) counterterms /ti(yu) and /t5(yu). With isospin 
conserved qI'"^ do not depend on K5, and furthermore are proportional to m^/(47r/)^, so 
these counterterms are small. Expressions for Fi and g'^ are given in the Appendix. 

The decay D* — > Dgn'^ is isospin violating, and the leading contribution occurs through 
7] — Ti^ mixing |^T|]. To first order in the isospin violation the decay is suppressed at tree 



level by the mixing angle 6 = (1.00 ± 0.05) x 10"^ M] 



1 _ (m„ - mrf) _ 2 ^ ^ 1_ 

Zl 2{ms-rh) 73 87.0' ^ ' 

Beyond tree level we have corrections to the rj—ir^ mixing angle parameterized by 6mix = 0.11 
p5| (Fig. 1^), loop corrections to the t] — mixing graph (Figs. Hb,c,d), as well as loop 
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a) 



D* 



D 



b) 
D* 



D 



c) 



D* / D* D* \ 



d) 



D 



D* 



I 1 

^ ' D 



FIG. 3. Nonzero vertex corrections for the decay D\ — > DgT:^ which involve tt^ — rj mixing. The 
cross denotes leading order mixing while the triangle denotes mixing at next to leading order. 

graphs with decay directly to vr" that occur in an isospin violating combination (Figs. ^,b). 
The contribution of Fig. ^ is again cancelled by the pseudo-goldstone boson wave function 
renormalization graph (Fig. ^d). Note that the decay D* DsTt^ cannot occur via a single 
virtual photon in the effective theory. In the quark model, decay to the spin and color singlet 
vr" can occur if the single photon is accompanied by at least two gluons (with suppression 
a/TT ^ 1/430 0]). We will neglect the possibility of such a single photon mediated transition 
here. Thus, 



1 



rrid 



2 (m<j — m) 



1 ~l~ ^mix ~l~ 



9^ 



^33 



log i^] 
mf 



+ -{-d^ + did2 + d^ 



2dido-2dl) +2Fi(m„di,rf2) -4Fi(mi,rfi,rfo)| 
2Fi{mi,di,d2) - A Fiijhi, di, do) 



(12) 



K^iog(|,; 



where for decay d^ = {6'''^ — 1)6, di = k-v + do, and d2 = — A + rfg. The tilde on the mass, 



rrii 



indicates that isospin violation is taken into account. Note that Y.i,b ^IbKb^ 
m\o. The function g^. depends on Ki, K5, and has both m\ and terms. 
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To describe D* D'-f, electromagnetic effects must be included, so the Lagrangian 
in Eq. (||) is gauged with a U{1) photon field -B'^. With octet and singlet charges, Q = 
diag(|, — |, — |) and Q' = | (for the c), the covariant derivative is defined as V^C, = 
d^^ + ieBij^[Q,C,] and V^H = d^H + ieB^^{Q'H — HQ) — V^H, where the vector and axial 
vector currents are now = \{i'^'D^^ + ^V^j,^'^) and = ^{^^V^^ — ^V^C,''). However, 
this procedure does not induce a coupling between D*, D and without additional pions. 
Gauge invariant contact terms should also be included, and it is one of these that gives rise 
to the D*D'~f coupling (and a D*D*'y coupling) 



Cf3 = ^TTHaH,a^''F^,Ql. 



(13) 
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Here (3 has mass dimension —1, = \{C^QC + CQ^^), and F^,^ = d^^B^, — d^B^. The terms 
which correct this Lagrangian at order uiq ~ 1/ rrtc have a similar form to those in Eq. (0) 

X X 
X X 



tuqA mqA 



e 



^ -g'Tri?,a^^/7,F,,, + ... . (14) 

The elhpses denote terms that do not contribute for processes without additional pions 
and/or can be eliminated using the equations of motion for H. For our purposes and 

in Eq. ([T^ ) are diagonal so only a = ai + a'l contributes. The finite part of as will 
be absorbed into the definition of (3. For mq = rric, the f3i term can be absorbed, and we 
absorb the part of the j32 term that contributes to D*D'y since D*D*'y only contributes in 
loops for us. Thus, (3 is defined to be the D*D'~f coupling at order l/rric. The last term in 
Eq. (JT^) is the contribution from the photon coupling to the c quark and has a coefficient 
which is fixed by heavy quark symmetry [|1|]. The ri^2 terms are similar to the ^2,3 terms 
in Eq. (0), and appear with /3 in the combination /? — (ri + T2)v ■ k/A"^. Here ti + T2 will 
have an infinite part necessary for the one-loop renormalization. Again it is not possible to 
isolate the finite part of the (ri + T2) contribution from that of f3, so the extraction at this 
order includes the renormalized 2 with v ■ k 0.137 GeV. 

The decays D*^ D°7, D*+ and ^ have decay rates F^, F^, and F^ 

given by 

K='i \^^a? = Zl, f/5 % + ^1 , (15) 
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where a ^ 1/137, k^ is the three momentum of the outgoing photon, and the wavefunction 
renormalization, Z^f, is given by Eq. (^. To predict the D* widths, Z'^^jZ'^ is kept to 
order and we take x 1/mc = l/mc. The vertex correction factor Z;^ has nonzero 
contributions from the graphs in Fig. ^ Note that the two one-loop graphs that contain 
a D^*^D*'K'y vertex (not shown) do not contribute |2y]. Furthermore, the graph in Fig. ^ 
has no contribution from the D*D*'~f coupling which arises from gauging the lowest order 
Lagrangian in Eq. (HI). Thus 



10 
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c) 



d) 
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D* D* S D* D 
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D 



FIG. 4. Nonzero vertex corrections for the decays D* D-y. 



1 + 



9^ KbQbb^ 



J log (-^) + + d,d2 + dl- 2d, do - 2 dl] 

mf ^ 



+2Fi{mi,di,d2) - 4 Fi(mi, rfi, c/q 



2 
3' 



(4vr/)2 2g, 



/i 

m,^ log ( — ^) + m; 
mf 



, 4^^ (A-,A;i)g- 
+^"(ai,tt5) , 



- log ( — )(rfo + + F2{mi, do, k ■ v) 



(16) 



where is the charge of meson vr*, is now the outgoing photon momentum, and the di 
are as above (again they differ depending on whether it is D* or one of D*^, D*^ that is 
decaying). The coefficients of the mj log {fi"^ / mf) terms agree with pO[. The new function 
F2 has mass dimension 1. It contains an analytic part proportional to 2do + v ■ k, and a 
non-analytic part which is a function of 6/mi and v ■ k/mi. q'^ contains the dependence of 
the rate on the (renormalized) counterterms Q;i(/i) and 0:5 (yu). Expressions for F2 and f?" are 
given in the Appendix. 

By examining Eqs. ([To|) , (|l3|), and (|l^) we can get an idea of the size of the various 
one-loop corrections to and F^*. With our power counting A ~ 5 ~ f ■ 



~ mq ~ ml so 



we can consider expanding in A/mj, 5/mi, and v ■ k/mi. Using the expressions from the 
Appendix gives 



Gi{mi,b) 



mj 



1 + T- + 



2 

mf 



Fi{mi,b,c) = 



1 - 



mi m. 
n{b 



c) I6{b^ + bc + c^) 



F2{mi, do, k ■ v) 



-Tcm^ 



mi 9mf 
3dl + Mok ■v + (k-v)"^ 
6mf 



+ 



(17) 



The leading terms in Gi and Fi are mq corrections to the rates. The second terms are order 
rrig/^ and ^Jm^jmc, and can be kept since they are unambiguously determined at the order 
we are working. The third and remaining terms in G\ and F\ are subleading in our power 



counting. The term —irmi in F2 is the formally enhanced contribution discovered in [|18 . 
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Note that there are no contributions to F2 proportional to 5 or k - v. The second term in F2 
in Eq. (|T7D has contributions of order m^^"^, Wmgk ■ v, and {k ■ vY/ Jm^ which again can be 



kept since they are unambiguously determined. 

The above power counting is sensible when is vhk or m^. We know that numerically 
~A~5~/i;-f, so for mj = the series in Eq. (|T^ are not sensible. In 
the term —nm.,^ in F2 was found to be important, so we want to keep corrections with 
dependence. Therefore, instead of expanding the non-analytic functions we choose to keep 
them in the non-analytic forms given in the Appendix. Numerically the one-loop corrections 
to and are very small; with g = 1 they are of order ~ 2%. For F^, 6mix is a 11% 
correction to the tree level result in Eq. (pT]) . Individually the terms proportional to g'^Fi 
and g"^ log {fi / rriq) in Eq. ([T3| ) are ~ 10% corrections for g = 1. However, the loops graphs 
with f] — mixing tend to cancel those without 77 — 7r° mixing leaving a ~ 2% correction. 
The one-loop corrections to F^ are larger, for instance the graph in Fig. ^ gives sizeable 
corrections that are not suppressed by g^. Corrections to the coefficient of the leading g'^/(3 
term obtained in |1^ range from ~ 3% for D* and ~ 20% for D*^ decay, to ~ 50% for the 
D*^. (The latter percentage is large because the only contribution for this decay come from 
a charged pion in the loop of Fig. ^d.) Corrections proportional to g"^ are only sizeable for 
D* Dg'-f where they are ~ 10% for g = 1. 



III. EXTRACTION OF g AND /? 

Using the calculation of the decay rates from the previous section, the couplings g and P 
can be extracted from a fit to the experimental data. Input parameters include rric =1.4 GeV 



27| , the meson masses from A = niD* — rriD = 0.142 GeV, 6 = rn^(*) — m^ji*) = 
0.100 GeV, and v ■ k which is determined from the masses. When isospin is assumed we 
use rriK = 0.4957 GeV and rrij^ = 0.1373 GeV. / is extracted from 7i~ decays. At tree level 
we use f = f^, = 0.131 GeV [0, while when loop contributions are included we use the 
one-loop relation between / and to get / = 0.120 GeV. The ratio of the decay rates 

F^ and F^ are fit to the experimental numbers 

13{D*^ D^-f)/B{D*^ DV) = 0.616 ± 0.076 [13] , 
B{D*+ D+-f)/B{D*+ D+7r°) = 0.055 ± 0.017 [0 , 

B{Dl ^ Ds7v^)/B{Dl Ds^) = 0.062 ± 0.029 [g , (18) 
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Order 


9 


/3(GeV-i) 




9 


/3(GeV-i) 




tree level 




= 3.6 


30. 








+Q' /rric + one-loop with 


0.23 


0.89 


4.3 


0.45 


2.8 


3.7 


+ chiral logs 


0.25 


0.78 


4.1 


0.56 


3.2 


1.4 


one- loop with nonzero A,(^,f • k, 


0.25 


0.86 


3.9 


0.83 


6.0 


2.5 


without analytic niq terms 














order ruq ~ l/rric with 


0.265 


0.85 


3.0 


0.756 


4.9 


3.9 


'^i = /^s = <5i = as = 















TABLE I. Solutions for g and /3 which minimize the \^ associated with a fit to the three ratios 



in Eq. (18). There are two solutions in the region of interest. 



where the errors combine both statistical and systematic. Using the masses m^j.o, m£,*+, 
rriD*, and mass sphttings m£)*o— m^io, m£)*+—mD+, mD'—Tnos from [|I^] gives the momentum 
ratios that appear in r"/r°: 

\kl\^ \klf |^3|3 



32.65 ±0.44, ^^ = 45.2 ±1.0, = 24.4 ± 1.5 . (19) 

|^l|3 |^2|3 |^3|3 ^ ' 

The errors here are clearly dominated by those in Eq. (|T8|) . Equating the numbers in Eq. (p!8D 
to the ratio of rates from Eqs. (|^) and (|TB|) gives a set of three nonlinear equations for g and 
/3 (where we ignore for the moment the unknown counterterms) . In general any pair of these 
equations will have several possible solutions. To find the best solution we take the error 
from Eq. ([T8|) and minimize the for the fit to the three measurements. We will restrict 
ourself to the interesting range of values, < g < 1 and < /3 < 6, discarding any solutions 
that lie outside this range. (The sign of g will not be determined here since it only appears 
quadratically in and F^J.) 

To test the consistency of the chiral expansion we will first check how the extraction of 
g and (3 differs at various orders. The results are given in Table |. At tree level only the 
ratio f3/g is determined, and the is rather large. We might next consider adding the 
contribution from the chiral loop corrections to D* —>■ D'j which go as y/m^. However, this 
does not lead to a consistent solution between the three data points unless (3 is negative. 
This signals the importance of the Q' /rric contribution in Eq. (P!3|) corresponding to a nonzero 
heavy quark magnetic moment. Adding this contribution gives the results in the second row 
of Table |I|, where there are now two solutions with similar in the region of interest. Adding 
the chiral logarithms, log {fi/rriq), at scale fi = 1 GeV gives the solutions in the third row. 
Taking nonzero 6, A, and v ■ k in the non-analytic functions Fi and F2 gives the solutions in 
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FIG. 5. Solution contours in the g-(3 plane for the situation in row 5 of Table |[ The solid, 
dashed, and dotted lines correspond to solution lines for the D*^, D*^, and D* decay rate ratios 
respectively. The stars correspond to the minimal solutions and the shaded regions correspond 
to the 68% confidence level of experimental error in the fit. The hatched region is excluded by the 



experimental limit T{D*+) < 0.13 MeV [|8|. 



the fourth row of Table |, where the value of g in the second solution has increased by ~ 50%. 
For these two solutions only the analytic mf dependence has been neglected. Finally, the 
solutions in row five include the analytic mf dependence with the counterterms set to zero 
(at fi = IGeV). The uncertainty associated with these counterterms will be investigated 
below. It is interesting to note that the extracted value of g in the second column of Table | 
changes very little with the addition of the various corrections. 

One can see more clearly how these solutions are determined by looking at Fig. |^. The 
central value for each ratio of decay rates in Eq. (|TBp gives a possible contour in the g-(3 
plane, as shown by the solid (D*^), dashed {D*'^), and dotted (-D*) lines. An exact solution 
for two of the ratios occurs at the intersection of two of these contour lines. However, a good 
solution for all three ratios requires a point that is close to all three lines. The solutions 
in the fifth row of Table | are indicated by stars in Fig. The size of the experimental 
uncertainties can be seen in the 68% confidence level ellipses which are shown as shaded 
regions in the figure (for two degrees of freedom they correspond to ^ Xmin + 2.3). These 
regions are centered on the solid line since the D*^ ratio has the smallest experimental error. 
The errors in Eq. ( pISD give the following one sigma errors on the two solutions 
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^7 = 0.2651:°?^ P = 0.85tfoGeV-\ ^ = 0.756l°:°^8 ^ = 4.9o+|7 CeV^ (20) 

Both solutions fit the first two ratios in Eq. ([T8|) , but do not do as well for the third. 
Minimizing the has biased against the third ratio as a result of its large experimental 
error. For this ratio the g = 0.265 and g = 0.76 solutions give values which are 4 and 
13 times too small respectively. For the first solution it is possible to improve the fit to 
the third ratio with reasonably sized counterterms. For instance, simply taking ai = 2 
gives B{D* Ds'K^)/B{D* Dsl) = 0.036. As we will see below, a large g solution with 

^ 1 is only possible if g increases to ~ 0.9 and /? increases to ~ 6.0 GeV^^ (c.f. Fig. ^. 

The experimental limit r{D*^) < 0.13 MeV translates into an upper bound on the 
value of g. Since B{D*^ D^l) is small, this bound is almost /? independent and to a 
good approximation is 

g < 0.52\/Vl + 3.01X- 1 a; = r(D*+) ""^'7(0.13 MeV) . (21) 

For the situation in row five of Table |I| this excludes the hatched region in Fig. ^. The limit 
on T{D*^) therefore eliminates the g ~ 0.76 solution at the two sigma level. Since this limit 
has not been confirmed by other groups it would be useful to have further experimental 
evidence that could exclude this solution. 

The central values in Eq. (^) have uncertainty associated with the parameter m^. Tak- 
ing = 1.4 ± 0.1 GeV gives 0.25 < g < 0.28 and 0.79 GeV"^ < (3 < 0.93 GeV"^ for 
the first solution, and 0.72 < 5^ < 0.80 and 4.6 GeV"^ < 13 < 5.3 GeV"^ for the second 
solution (in both cases the changes very little). There is also ambiguity in the solu- 
tion in Eq. ( pOD due to the choice of scale (ie., the value of the counterterms ai, cts, 
ki and K5). Increasing /i to 1.3 GeV gives solutions {g = 0.28, /3 = 0.91GeV~^,x^ = 1-4) 
and {g = 0.78,/? = 5.0GeV^^,x^ = 4.1), while decreasing /i to 0.7 GeV gives solutions 
{g = 0.25, /3 = 0.83GeV-\x2 = 3.7) and {g = 0.72, p = 4.7GeV~\x^ = 3.1). Note that 
the of the second solution remains large, while the of the first solution is reduced 
significantly by an increased scale. 

Another method of testing the effect of the unknown counterterms ai, a^, ki and K5 
is to take their values ai fi = 1 GeV to be randomly distributed within some reasonable 
range of values. We take —1 < iii,K,^ < 1 and —2 < ai,a^ < 2, with the motivation 
that the counterterms change the tree level value of Z!^ and Z!^ by less than 30%, and give 
corrections that are not much bigger than those from the one-loop graphs. Near each of the 
two solutions 5000 values of g and (3 were then generated by minimizing the x^- This gives 
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2 4 6 8 10 

p (GeV"') 

FIG. 6. Effect of the order niq counterterms k^, ai, and as) on the solutions in Eq. (| 
The counterterms are taken to be randomly distributed with —1 < ki,k^ < 1, —2 < 01,05 < 2. 
For each set of counterterms g and (3 were determined at the new minimal x^- 5000 sets were 
generated near each of the two solutions. 

the distributions in Fig. The solution with g = 0.265 and (3 = 0.85 GeV^^ has fairly small 
uncertainty from the counterterms. The g = 0.76, (3 = 4.9 GeV~^ solution has much larger 
uncertainty because the corresponding contour lines in Fig. § are almost parallel. For this 
solution the upper bounds are determined by the limits of a few MeV on the D* widths. 
From this analysis we estimate the theoretical uncertainty of the solutions in Eq. ( ]20|) to be 
roughly 



g = 0.2651°;°^ (3 = 0.851°;? GeV"^ 



g = 0.76l°;2 (3 = 4.9l^:° GeV"^ 



(22) 



at this order in chiral perturbation theory. The errors on g and f3 are positively correlated 
since the values of g and (3 are constrained in one direction by the small error on the 
rate ratio in Eq. (P^). 

From Eq. (|l|) and Fig. |, we see that if the error in B{Dl D^tx^) /B{D* 1^,7) 
can be decreased by a factor of two, in conjunction with a limit of r(D*+) ^ 0.6 MeV then 
this could provide strong evidence that the g = 0.76 solution is excluded. On the other 
hand if the central values of the second and third ratios in Eq. ([T8| ) decrease, then a width 
measurement or stronger limit on r{D*^) will be needed to distinguish the two solutions. 

Using the extracted values of g and f3 gives the widths shown in Table |I[ The couplings 
were extracted at one-loop and order ~ ^/rric, so the predictions for the D* widths are 
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Predicted widths in keV 






D* 

^ s 


B*+ B*° B* 


g = (J. zoo, p = (J.ooGeV 


18 


26 


0.06 


~ U.Oo ~ U.03 ~ 0.04 


uncertainty from experiment 


16 - 24 


23 - 35 


0.01 - 0.13 




uncertainty from counterterms 


16 - 27 


22 - 39 


0.04 - 0.13 




5 = 0.76, /3 = 4.9GeV-i 


323 


448 


103 


~ 2.1 ~ 2.0 ~ 1.6 


uncertainty from experiment 


285 - 367 


396 - 508 


83 - 128 




uncertainty from counterterms 


215 - 1318 


281 - 1157 


53 - 1078 





TABLE II. Widths in keV for the D* and B* mesons. The experimental and counterterm 
ranges are determined by the extremal values of g and (3 in Eqs. ( pO| ) and (p^). For g = 0.265 the 
D* width is small due to a delicate cancellation in /is as explained in the text. The uncertainty in 
the B* widths is large due to unknown l/rric^b corrections. 

made at this order. The experimental uncertainty in the D* widths is estimated by setting 
g and (3 to the extremal values in Eq. (pUD, which gives the range shown in the second and 
fourth rows of the Table. The uncertainty from the unknown counterterms in the third and 
fifth rows is estimated in the same way using the uncertainties from Eq. (^). Note that for 
the g = 0.265 solution the D* width is small due to a delicate cancellation in fi^ resulting 
from setting Z^^ x l/rric = l/rric. Keeping Z^f/rric to order gives a D* width of 0.28 keV 
with a range of 0.1 — 0.4 keV for both the experimental and the counterterm uncertainties. 

Making use of HQS allows us to predict the width of the B* mesons from their dominant 
mode B* B'-f. Eq. ([T5|) gives the rate for B* B^ with Q' = —1/3 and rric rrib. Since 
the couplings /?i and /?2 are unknown these rates can not be determined at order l/m^b, but 



we can include the order rUg corrections. The B meson masses are taken from |0 and we 
use rrib = 4.8 GeV ||2^. We set 6 = 0.047 GeV and A = k ■ v = 0, but since the contribution 
Q'/rrib in Eq. (|r^) is numerically important it is kept in our estimate. For comparison the 
widths obtained with the g = 0.76 and f3 = 4.9 GeV~^ solution are also shown. 

As a final comment, we note that heavy meson chiral perturbation theory can also be 
used to examine excited D^*'^ mesons, such as the p-wave states, Dq, Dl, Di, and D2 



29| , pT|| . To do so, explicit fields for these particles may be added to the Lagrangian giving 



a new effective theory. For interactions without external excited mesons (such as the ones 
considered here) these new particles can then contribute as virtual particles. However, since 
we have not included these heavier particles they are assumed to be 'integrated out', whereby 
such contributions are absorbed into the definitions of our couphngs. 
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IV. CONCLUSION 



For the and D*, the decays D* Dtt and D* D'y are well described by 

heavy meson chiral perturbation theory. Using the recent measurement of B{D*^ D'^l) 
|l^, the ratios of the and Dti^ branching fractions were used to extract the D*D7t and 



D*Dj couplings g and /5. Two solutions were found 

^ = 0.2651:°^!:°^ P = 0.85titiGeV-' g = 0.76t-',lti P = A.9titfGeV-K (23) 

The first error here is the one sigma error associated with a minimized fit to the three 
experimental branching fraction ratios (see Fig. ^). The second error is our estimate of the 
uncertainty in the extraction due to four unknown counterterms ai, 0:5, ki and ^5 that arise 
at order niq (see Fig. H). 

It is possible that the uncertainty from these counterterms can be reduced by determining 
them from other processes. For these corrections to contribute at low enough order in the 
chiral expansion we need processes with outgoing photons or pseudo-goldstone bosons, such 
as semileptonic D decays to K, rj, or vr. Here there are also SU(3) corrections to the left 
handed current which involve an unknown parameter 7]q [0. Information on Ki and k[ can 
be determined from the pole part of the Dg — > Kiui form factor 0. In a similar manner 
Ds rjiui can constrain ki and K5, and a comparison of the form factors for — > K^iui 
and Ds rjiui gives information on k[ and K5. These investigations were beyond the scope 
of this paper. In principle, information about the constants ai, and 0:5 could be obtained 
from a measurement of B ^ 'yii^e- The CLEO experimental bound on B ^ Cvi {i = e, fi) 
Pm is roughly two orders of magnitude above the theoretical prediction, but due to the 



helicity suppression for B — > iu^ the branching ratio for B — > 'yiui may be up to an order of 



magnitude bigger |31 



Another possible approach would be to use large Nc scaling for the counterterms in 6Cg 
and 6Cf3. Terms that have two chiral traces are suppressed by a power of Nc compared to 
those with only one trace. In the large Nc limit the counterterms ki and ai would dominate, 
and ^5 and could be neglected, thus reducing the theoretical uncertainty. 

The smaller solution for g in Eq. ( p3D is fairly insensitive to the addition of the one-loop 
corrections (see Table 1). However, corrections at order rriq ~ ^/rric, including the heavy 
meson mass splittings, were important in determining the solution with larger g. The limit 
r(D*+) < 0.13 MeV [^] gives an upper bound on the coupling g (see Eq. (pl|) and Fig. 
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and eliminates the g = 0.76, (3 = 4.9 GeV~^ solution. Experimental confirmation of this limit 
is therefore desirable. Note that the largest experimental uncertainty in our extraction comes 
from the measurement of B{D* DsTt^), and dominates the theoretical uncertainty due to 
decay via single photon exchange. A better measurement of B{D* Ds'k^)/B{D* Dg'y) 
along with a limit r(D*+) ^ 0.6 MeV could provide further evidence that the g = 0.76 
solution is excluded. However, if the central values of the second and third ratios in Eq. (|TBp 
decrease then a width measurement or stronger limit on F will be needed to distinguish 
the two solutions. An improved measurement of B{D* DgTi^) may also give valuable 
information on the unknown couplings ki, k^, ai, and a^. 

The extraction of g has important consequences for other physical quantities [2-11]. For 
example ||32|, for the B niui form factors with < 2m.^, analyticity bounds combined 
with chiral perturbation theory give g fs ^ 50 MeV ||33|. The solution g = 0.265 gives 
fs ~ 190 MeV for the B decay constant. However, for g = 0.76 we have fs ^ 66 MeV, 



which is roughly a factor of three smaller than lattice QCD values, ~ 160 — 205 p4 
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APPENDIX: ONE LOOP CORRECTION FORMULAE 



In this appendix we give explicit formulas for the functions Gi, Fi, F2, p^, and g?^ that 
occur in our one loop correction formulae in Eqs. (|^), (0), ([T3|), and (p^Bf). In doing this 
type of one-loop calculation an important integral is 



2e 



ib 



^ + ln^ -f 2-2F - 



(27r)4-2^ (g2 -rn"^ ^is)2{q-v-b + le) {Any 
where l/e = 1/e — 7 + log (47r). F is needed for both positive and negative b, so 



(Al) 



X 



^ — tan ^ 



\x\ < 1 
\x\ > 1 



(A2) 
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For 6 > the function F was derived in ||35|Q and agrees with the above formula^. For 
X = h/m < —1 the logarithm in Eq. ( |A2|) has an imaginary part. This corresponds to 
the physical intermediate state where a heavy meson of mass mu produces particles of 
mass niH + h and m. For the calculation here the imaginary part only contributes from 
F(m,r/('^o A))) and was found to always be numerically insignificant. Note that the real 
part of xF{l/x) is continuous everywhere, and differentiable everywhere except x = — 1. 
Also F(l) = F(-l) = 0. 

Eq. d^) contains the function 
5 



Gi{a,b) = 



{Ab' -^a')F{a/b) + ^{a' -b')^F'{a/b) 



(A3) 



where is the analytic contribution. In the limit A — Eq. (j^) gives Zd = Zd* in 
agreement with HQS. To obtain HQS in the finite part of the dimensionally regularized 
calculation of the graphs in Fig. ^ it was necessary to continue the D* fields to d = 4 — 2e 
dimensions (so the D* polarization vector ea = (1 — |)ea where J^^a^"" = ~3). 
In Eqs. (0) and ([T3|) we have the functions 



b{a^ - b^)F{a/b) - c{a^ - ^)F{a/c) 



9l 



(4vr/)2 
1 



2 



(A4) 



rrid 



(47r/)2 2(m, - m) 



+ 



We have ignored isospin violating counterterm corrections in f)^'^ and work to leading order 
in the isospin violation for g^. In deriving Eq. ( |A^ ) use has been made of = 23^711^ = 



2Bomd = 2Borh, m\ — m^/2 = BqUIs, and m'^± — m^o = 
In Eq. (|16D we have the function Fi and the function 

2a2 



md)Bo. 



F2(a,b, c) 



-2 b 



-2 b 



c 

2a2 

c 



Fit) 

a/b t'^ 



1 

4x2 



F{x) F{x) 



2x2 



4(x2 - 11 



x=a/ (fe+c) 



x=a/b 



(A5) 



Assuming isospin to be conserved the counterterm contributions in Eq. (|T6D are 



^Eq. ( |A2D for F disagrees with ^ for a; < 0. Their F{l/x) is even under x — > —x making Eq. (Al) 
discontinuous at A = 0. Furthermore, their F has no imaginary part corresponding to the physical 
intermediate state. 
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ml ai _ {m\ - ml) 
2(47r/)2 3Q,„(47r/)2 ' 
(2m|^ - ml) Oil {m\ ~ ml) 
2 (47r/)2 + (47r/)2 
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